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Rotational Correlation Functions of Single Molecules 
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Single molecule rotational correlation functions are analyzed for several reorientation geometries. 
Even for the simplest model of isotropic rotational diffusion our findings predict non-exponential 
correlation functions to be observed by polarization sensitive single molecule fluorescence microscopy. 
This may have a deep impact on interpreting the results of molecular reorientation measurements 
in heterogeneous environments. 
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Reorientation of small molecules or segments of macro- 
molecules undergoing conformational changes are ele- 
mentary processes and often of crucial importance for 
the properties of materials. In most of the experimen- 
tal techniques probing rotational dynamics, ensembles of 
molecules have been monitored 

HSU- Actually, many 
chemical and biological systems are heterogeneous, a fact 
that e.g. is reflected in their dynamical properties. In the 
context of the non-exponential decay of bulk correlations 
in supercooled liquids this has been termed dynamic het- 
erogeneity. Proteins, polymers, colloids and coarsening 
systems share this feature of non-exponential relaxation. 
In the past, bulk techniques have been developed that 
allow a characterization of heterogeneous rotational dy- 
namics in complex systems in terms of a distribution of 
reorientational rates fluctuating in timeHQ. There arc, 
however, hints that even without (bulk) heterogeneity 
rotational correlations functions (RCFs) show deviations 
from exponential decay on a microscopic scale 6] . Such 
a behavior has also been found in numerical calculations 
on spin glasses Q, • 

In order to get a deeper insight into the nature of the 
rotational dynamics and its relation to structural prop- 
erties of heterogeneous materials, it is highly desirable to 
conduct experiments providing local (as opposed to aver- 
aged bulk) information. It has been shown by numerous 
examples that a single molecule is an exquisitely sensitive 
probe of dynamical processes in its local environment 0, 
llCj . Along these lines, temporal fluctuations of excited 
state lifetimes [TH. excitation^^ an( A emission spectra[l3j 
and electronic coupling strengths[lj] have been studied 
by single molecule spectroscopy and analyzed in terms 
of fluctuating environments. Here we focus on the ro- 
tational dynamics of single fluorophores investi gate d b y 
polarization resolved microscopy [15|, Il6l Il7l Il8l IliA |20|. 
In order to extract quantitative and reliable informa- 
tion from such measurements, a detailed analysis of local 
RCFs is of utmost importance. It is the purpose of the 
present letter to provide such an analysis for the simple 
model of rotational diffusion of single fluorophores. 



Fluorophores preferentially absorb photons whose elec- 
tric field vectors are aligned parallel to the transition 
dipoles of the molecule 21] . These dipoles have a well 
defined orientation with respect to the molecular axes. 
Similarly, emission also occurs with light polarized along 
a fixed axis. The experimental situation we will ana- 
lyze in the remainder of this paper is depicted in Fig. 1 
together with the definition of the axes system and the 
angles used. The orthogonal polarization resolved emis- 
sion signals I p and I s depend on both, absorption and 
emission efficiency. The contributions of I p and I s to the 
signals collected in fluorescence correlation spectroscopy 
have been analyzed in detail [12, H3- 

A customary measure of dipole orientation can be ob- 
tained by calculating the reduced linear dichroism d from 
I p and I s via 



(1) 



This frequently used quantity[H E H El 

fluctu- 
ates in the course of time due to molecular rotational dy- 
namics. The normalization ensures that d depends solely 
on the orientation of the emission dipole and not on the 
absorption efficiency. Hence we represent the transition 
dipoles not by vectors but by their orientation using an- 
gles 9, 4>. We want to point out that the analogy of d 
to what is often called polarization[2l| could lead to the 
wrong assumption to use the anisotropy values instead. 
However, in contrast to fluorescence depolarization meth- 
ods this would not make any sense here. Subsequently 
the auto-correlation function of d(t) is calculated, 



C d {t) = (d(t + T)d(r)) 



(2) 
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In the present letter we will show that Cdit) is a non- 
trivial correlation function which nevertheless gives valu- 
able information about the time scale and the geometric 
aspects of molecular reorientations. It should be noted 
that in principle the fraction I p /I s leads to the azimuth 
<fi, however, an extremely well S/N ratio would be needed 
for reliable values. 

Well known bulk techniques allowing to monitor molec- 
ular rotational dynamics are provided by dynamic light 
scattering 0, dielectric spectroscopy Q or nuclear mag- 
netic resonance (NMR) j3(. Due to the different ranks I 
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FIG. 1: Schematic experimental setup: The polarization of 
light emitted by a single fluorophore depends on the orien- 
tation of the transition dipoles with respect to the detection 
system. The two orthogonally polarized components of the 
signal are separated by a polarizing cube beam splitter. 



of the respective interaction tensors distinct RCFs 
d(t) = (A(cos(0(t))fl(cos(0(O))) 



(3) 



of the relevant Legendre polynomial are obtained. Here, 
9 is defined in the same way as in Fig.l. For instance, 
the interaction of a permanent electric dipole moment 
with an applied electric field transforms in the same way 
as the first order Legendre polynomial Pi (cos 9) — cos 9 
and therefore C\ (t) is observed in dielectric spectroscopy. 
NMR and light scattering are I = 2 techniques and yield 
C2(t). Below it will be shown that Cd(t), apart from 
being a local correlation function, cannot be related to 
one of the Ci(t) in a simple manner. 

The theoretical description of reorientational dynamics 
in molecular liquids is a complicated many particle prob- 
lem. For a qualitative understanding of the behavior of 
rotational correlation functions, however, a treatment in 
terms of rotational Brownian motion is usually sufficient. 
In the present letter we restrict ourselves to the simple 
model of anisotropic rotational diffusion, thus neglecting 
inertial effects completely. This approximation is well 
suited for the systems we have in mind, namely biologi- 
cal systems or supercooled liquids and glasses, where the 
time scale of the relevant relaxational modes is slow com- 
pared to typical vibrational frequencies. 

In all following considerations we assume the transition 
dipole matrix elements to be stationary quantities. In 
this case we have Cd(t) — (d(t)d(0)), with 



d{9,4>) = 



cos 2 9 — sin 2 9 cos 2 i 
cos 2 9 + sin 2 9 cos 2 i 



(4) 



Note that d(9, 4>) cannot be expressed in terms of the 
spherical harmonics Y2 m {9, </>) and this is why Cd(t) can- 
not be related to the corresponding RCF C^)- This is 
to be contrasted to the optical anisotropy, the correlation 
function of which is directly proportional to (t) . 

In order to keep the treatment general, in the following 
we treat auto-correlation functions of the form 



c x (t) = (x(Q PL (t))x(n PL (o))). 



(5) 



Here, X denotes some arbitrary function of the molecular 
orientation Qpl of the principal axis system (P) of the 



relevant interaction tensor in the laboratory fixed frame 
(L) . In the present context we are primarily interested in 
the case of fluorescence, X = d, in which case the z-axis of 
the P-system coincides with the direction of the transi- 
tion dipole. However, other choices like e.g. X = e iqr rel- 
evant for incoherent neutron scattering can be analyzed 
in an identical way. If the molecular orientation O(t), 
given in terms of Eulerian angles |25|. is modelled as a 
stochastic process, Cx(t) can be written in the form[26l| : 



C x (t) = ^ J dn J dfloX(Q)X(Q )P(Q, t\Q ) 



(6) 



where -P(0, i|0o) denotes the conditional probability to 
find f2 at time t given Qq at time t — 0. Throughout this 
letter we restrict ourselves to the model of anisotropic 
rotational diffusion of symmetric top molecules, in which 
case one has 



e -[i((+l).D x +m 2 (-D,--D*)]t 



(7) 



Here, the 
elements^ 



are Wigner rotation matrix 
The rotational diffusion coefficients D y 
and D x are equal, but different from D z . The limit of 
isotropic rotational diffusion is recovered for D x —D z . It 
is important to point out that O in eq.Q denotes the 
orientation of the diffusion tensor (D) in the L-system, 

In order to proceed in the calculation of Cx(t), 
we expand the quantity X(£Ipl) in the expres- 
sion for the correlation function, eq.jtjj), in terms of 



D 



(0 



with 



X((l PL (t)) = Y J l, m n X hmnD^L^lpL{t)) 

Xi- mn = &£fd£lD$Z(Sl)X(Si). Next, D%l(n PL (t)) 
is expressed in terms of the relevant Od.l(£) via 
= fipc + £l DL {t) using D§n(Sl P L,{t)) = 
J2n Dnl(^PD)D i l ll{n DL (t)). Here, Q PD = (a, 9 pd ,tt - 
4>pd) is a shape and symmetry dependent molecular 
quantity. For most relevant cases X(VLpp>) has axial sym- 
metry, which allows to replace the Dq^ by spherical har- 
monics Yim- Performing the calculation using eqns.© 
and |0 yields: 



C x (t) = J2 A i C i( t )- 



With A, = i£ m ^d<\> JZd9 sin 9X (9, <j>)Y lm ( 
and 



(8) 



Q(t)=E|€w 



-[;(i+l)_D x +m 2 (_D,-_D x )]t 



(9) 



the RCFs eq. (JHJ are calculated numerically Typically, it 
is sufficient to use I values up to I = 20. Thus, Cd{t) is 
given by a weighted sum of a large number of RCFs. This 
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FIG. 2: Rotational correlation function Cd(t) for isotropic 
diffusion (solid line). For comparison the exponential corre- 
lation function I — 2 (ea llOl A2 = X,Ai^2 = 0) is plotted as 
dotted line 

fact has not been noticed in the literature before to the 
best of our knowledge. As already mentioned above, for 
the optical anisotropy one has (r(i)r(O)) = C 2 (i) which 
decays as a superposition of at most three exponentials. 

For the case of isotropic rotational diffusion, i.e. D = 
D x = D y = D z , cq.@ simplifies to 

cx(t) = j2 A ^ l{l+1)Dt ( 10 ) 
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FIG. 3: (a) Stretching exponent /3 depending on the 
anisotropy S — D z /D x (D x = D y ). Qpd denotes the angle 
between the z-axis of the diffusion tensor and the transition 
dipole ft. (b) Correlation times normalized by to — t(D z || ft). 
t and (3 were obtained by fitting C(t) with exp( — (t/r) ,J ) 



For symmetry reasons, the odd components vanish. A 
similar expression has been reported recently [24]. In 
table [I] the pre-factors are listed up to I = 20. The re- 
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TABLE I: Pre-factors Ai from eq.dJJ for isotropic rotational 
diffusion 

suiting correlation function significantly deviates from 
an exponential decay as shown in Fig. 2. From fitting 
the computed data with a stretched exponential function 
f(t) = e -(*/ T <=)' 3 we obtain r c = 0.87/6Z? and (3 = 0.871. 
Therefore, even for this simple model we find intrinsic 
non-exponential relaxation. The optical anisotropy de- 
cays exponentially, (r(i)r(O)) —e^ eDt . 

The model of stochastic isotropic reorientations 
presents a drastic approximation as most molecules do 
not exhibit a spherical shape. Thus, in many cases of 
practical interest one is confronted with asymmetric top 
molecules. As noted above, however, in the present let- 
ter we restrict ourselves to the case of symmetric tops for 
simplicity. This means that we consider an ellipsoidal dif- 
fusion tensor with three components D x — D y and D z , 
see Fig. For standard I = 1,2 methods the RCFs 



resulting for this model consist of two and three expo- 
nentials, respectively, with time constants depending on 
the degree of anisotropy 5 — D z /D x . It is evident from 
eq. @ that for D x ^ D z the degree of non-exponentiality 
is more pronounced due to the fact that for each I now 
one has I + 1 exponentially decaying functions with dif- 
ferent weights. However, the orientation of the transition 
dipoles with respect to D z determines the degree of this 
additional non-exponentiality. Only for D z \\ fl the same 
stretching exponent as for isotropic reorientation is ob- 
tained. 

An entirely different approach to the calculation of 
RCFs via rotational random walk simulations has proved 
to be successful in the context of supercooled liauids[27|. 
Here, molecular reorientation has been modelled by rota- 
tional jump processes. Starting from an arbitrary orien- 
tation, after a certain waiting time drawn from Poisson 
statistics the next orientation is chosen at random with 
the only restriction that starting and end position differ 
by an angle 7. In the limit of 7 — > isotropic rota- 
tional diffusion is obtained with results identical to those 
obtained with the analytic approach. In Fig.0|the influ- 
ence of the jump angle 7 on the correlation function Cd(t) 
defined by eqns.QJ and © is plotted. With increasing 7 
the correlation decays become more exponential, which 
essentially originates from an convergence of the different 
I contributions to the RCF. 

So far, we have assumed that the principal axes sys- 
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FIG. 4: Influence of the reorientation mechanism on the 
stretching (o) and the time scale (•) of the correlation func- 
tion, calculated for rotational jump processes with varying 
jump angles 7. For comparison the correlation times were 
normalized by the pure / = 2 correlation time T2. 

tem (P) relevant for experimental observables has a fixed 
orientation flpo relative to the 'diffusion tensor' sys- 
tem (D), and only the orientation of the latter with re- 
spect to the laboratory system (L) is time dependent: 
f2(i) = f2rji(t). In principle, one can treat the case 
of internal rotations by assuming a composite Markov 
process |26( {Qpr)(t), £lDL(t)} and using an appropriate 
master equation for the composite process. If the inter- 
nal motion is independent of the tumbling motion of the 
whole molecule one can factorize the corresponding prob- 
ability functions and averages. This case is of particular 
relevance for the conformational motions of polymers or 
proteins in solution and has been treated extensively in 
the literature to which we refer [2^. 

In the present analysis we have neglected the influence 
of the numerical aperture (NA) on the optical path. Us- 



ing very high NA objectives (NA > 1) out-of-plane con- 
tributions of the transition (emission) dip ole to the flu- 
orescence signal have to be considered [la fl8| . However, 
it was pointed out that even for NA= 1.25 the dichroism 
signal is only slightly influenced by this effect [18]. Addi- 
tionally, if one were to measure rotational dynamics over 
an extended temperature range, it is experimentally more 
convenient to use low NA (< 1) objectives. Here still an 
influence of the NA is expected, which, however, becomes 
smaller the smaller the NA is. Nevertheless, the accurate 
incorporation of the NA into our theoretical description 
has still to be worked out. We also would like to men- 
tion that methods have been developed which utilize the 
longitudinal field component of high NA objectives to ex- 
tract the full three dimensional orientation of transition 
dipoles 29, 30]. Although quite powerful, such measure- 
ments require a very good S /N ratio for subsequent data 
analysis and are limited with regard to time resolution. 

While in the past only bulk experimental techniques 
were available for studying molecular dynamics, recent 
experimental progress has enabled to study dynamical 
processes at the single molecule level. In particular, flu- 
orescence microscopy has allowed to track molecular ori- 
entation in time. Here the merits of single molecule mi- 
croscopy emerge in probing rotational dynamics in het- 
erogeneous environments in a direct way. We found that 
even for the model of isotropic rotational diffusion the 
obtained correlation function decays non-exponentially. 
Motional anisotropy increases the deviation from expo- 
nentiality. If the RCFs calculated for the simple model 
of anisotropic rotational diffusion are parametrized by a 
stretched exponential, we find stretching parameters (3 > 
0.85 for physically reasonable values of the anisotropy 8 
between 0.3 and 3. From this we conclude that only 
(3 values smaller than roughly 0.85 can be taken as 
an indication for intrinsically non-exponential dynam- 
ics. Our findings suggest a careful interpretation of single 
molecule rotational correlation functions as obtained by 
polarization resolved microscopy. 
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